Instanton for the Kraichnan passive scalar problem 
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We consider high-order correlation functions of the passive scalar in the Kraichnan model. Using 
the instanton formalism we find the scaling exponents of the structure functions Sn for n ^ 1 
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INTRODUCTION 

Anomalous scaling is probably the central problem of 
the theory of turbulence. In 1941 Kolmogorov formulated 
his famous theory of developed turbulence ||l| , where the 
scaling behavior of different correlation functions of the 
turbulent velocity was predicted. Experimentally one ob- 
serves deviations from the scaling exponents, proposed by 
Kolmogorov It is recognized that the deviations 

are related to rare strong fluctuations making the main 
contribution into the correlation functions [g|-|7|- This 
phenomenon, which is usually called intermittency, is the 
most striking peculiarity of developed turbulence. 

One of the classical objects in the theory of turbulence 
is a passive scalar advected by a fluid. The role of the 
passive scalar can be played by temperature or by the 
density of pollutants. Correlation functions of the scalar 
in a turbulent flow possess a scaling behavior that was 
established by Obukhov and Corrsin |^,^ in the frame 
of the theory analogous to that of Kolmogorov . Inter- 
mittency enforces deviations from the Obukhov-Corrsin 
exponents that appear to be even stronger than the devi- 
ations from the Kolmogorov exponents for the correlation 
functions of the velocity ||l0|-[l3|| . 

Unfortunately, a consistent theory of turbulence de- 
scribing anomalous scaling has not been constructed 
yet. This accounts for the difhculties associated with 
the strong coupling inherent to developed turbulence. 
This is the reason for attempts to examine the inter- 
mittency phenomenon in the framework of different sim- 
plified models. The most popular model used for this 
pur pose is Kraichnan's model of passive scalar advection 
|l4| , where the advecting velocity is believed to be short 
correlated in time and have Gaussian distribution. That 
allows one to examine the statistics of the passive scalar 
in more detail. 

The scalar in the Kraichnan model exhibits strong in- 
termittency even if it is absent in the advecting velocity 
field. This was proved both theoretically [T5|-|24| and nu- 
merically p3|-p5[|. In the theoretical works the equation 
for the n-point correlation function was solved assum- 



ing that different parameters, such as ^2, 2 — ^2, or d"^, 
are small (recall, that ^2 is the exponent of the second- 
order correlation function of the passive scalar and d is 
the dimensionality of space). The order of the correla- 
tion functions that can be examined in the framework of 
the methods of the noted papers is bounded from above, 
which does not allow one to imagine the whole depen- 
dence of Cn on n. For that it would be enough to get 
the asymptotic behavior of at n ^ 1. There have 
been several attempts to find the scaling of the correla- 
tion functions for larger n. In the the work by Kraichnan 
p6| a closure was assumed enabling one to find C„ for 
any n. An alternative scheme was proposed in ]27[ |. An 
attempt to solve the problem at large n was made in [ p8[ , 
where an n-independent asymptotic behavior was found. 

In the present work we develop a technique based on 
the path-integral representation of the dynamical corre- 
lation functions of classical fields ||2^-^^. We use an 
idea, formulated in ||3^, that is related to the possibil- 
ity of exploiting the saddle-point approximation in the 
path integral at large n. The saddle-point conditions are 
integro-differential equations describing an object that, 
in analogy to the quantum field theory, we call an in- 
stanton. The instantonic method was already success- 
fully used in some contexts. Results concerning Burgers 
turbulence, conventional Navier-Stokes turbulence, and 
modifications of the Kraichnan model were obtained with 
the help of this method in the Refs. The formal- 

ism presented in this paper enables one to find correla- 
tion functions of the passive scalar for arbitrary n 3> 1 
provided dC^2 3> 1. 

The paper is organized as follows. In Sec. | we for- 
mulate the Kraichnan model, introduce notation, and 
write down the standard path integral representation for 
the correlation functions. This basic representation turns 
out to be unsuitable for the saddle-point approximation; 
therefore, we reformulate the problem in Sec. |^ Pass- 
ing to new variables that are Lagrangian separations, we 
get a path integral that already admits the use of the 
saddle-point approximation. In Sec. [II we consider the 



instantonic equations for the case of the structure func- 
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tions. We solve these equations in the hmit rfC2 ^ 1, 
which enables us to find the anomalous scaling and esti- 
mate the n dependence of Sn ■ T he main results of the 
work are presented in Sec. Ill C and discussed in Con- 
clusion. Details of calculations are given in Appendixes. 



I. KRAICHNAN MODEL 

Advection of a passive scalar by a velocity field v is 
described by the equation 



(1.1) 



where k is the diffusion coefficient and is the source 
of the passive scalar (say, if corresponds to fluctu- 
ations of temperature, then cf) represents the power of 
heaters). In a turbulent flow, d is a random function of 
time and space coordinates. The source (j) is also assumed 
to be a random function. Then passive scalar correla- 
tion functions are determined by the statistics of v and 
<j). Usually, one is interested in simultaneous correlation 
functions F„ = {0{ri) ■ ■ ■ 9{rn)) since a large-scale veloc- 
ity destroys temporal correlations in the Eulerian frame, 
whereas simultaneous objects are not influenced by it. 

It is convenient to examine the anomalous scaling in 
terms of the structure functions 



S^{r) = {\e{r/2)-e{-r/2T) 



(1.2) 



One expects a universal behavior of the structure func- 
tions in the convective interval of scales ^ r <C i, 
where ra is the scale where the diffusivity becomes rele- 
vant and L is the correlation length of the scalar source 
(j). Namely, one observes a scaling dependence on r: 



Sn{r) cx / 



(1.3) 



In the frame of the Obukhov theory Cn = {n/2)C,2- 
Therefore, the differences {n/2)C,2 — Cm which are usually 
called anomalous exponents, characterize the anomalous 
scaling. One can write an estimate 



5„(r)^ A„ [^2(r)]"/' 



(n/2)C2-Cn 



(1.4) 



where An is an n-dependent factor. Note that Eq. 
(1.4) implies that the structure functions in the con- 
vective interval do not depend on the diffusion length 
Td- The intermittency leads to the conclusion that val- 
ues of the structure functions should be much larger 
than their naive Obukhov estimations 0. Therefore 
(n/2)C2 — Cn > Eind we conclude that these are the 
anomalous exponents that reflect the intermittency. 



A. Formulation of the problem 

In the Kraichnan model both v and cf) are assumed to 
be independent random functions, 5 correlated in time 



and described by Gaussian statistics homogeneous in 
space. Therefore, statistical properties of the fields are 
entirely characterized by the pair correlation functions 



(0(ii, ri)0(t2, r2)> = x{ri2)5{ti - ^2) , x(0) 
("0(^1, ri)u^(t2, r2)) = Vap{ri ~ r2)5{ti - 



P2 



Here, xi'"') is a smooth function decaying on the scale L 
which is the pumping length. The constant P2 has the 
meaning of the pumping rate of 9^. The tensor Vapir) 
has a characteristic scale L„, which has the meaning of 
the pumping length of the velocity. We will assume that 
3> L. Since r <C L in the convective interval, we will 
need Vap only dX r <^ L^, where one can write 



(1.5) 



The quantity Vq is an r-independent constant that is 
the main contribution to the velocity correlation func- 
tion on scales less than the velocity pumping length L^. 
Nevertheless, besides Vo, we should also keep a small r- 
dependent correction /C, since Vq corresponds to advec- 
tion homogeneous in space and therefore does not con- 
tribute to simultaneous correlation functions of 6. 

The velocity correlation function is assumed to possess 
some scaling properties. Namely, /C(r) oc r^"''', where the 
exponent 7 characterizes the roughening degree of the ve- 
locity field. The field is smooth in space at 7 = and 
is extremely irregular at 7 = 2. We will treat an arbi- 
trary 7 satisfying the inequality < 7 < 2. The tensorial 
structure of ICafs is determined by the incompressibility 
condition diw = 0, implied in the Kraichnan model 
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^_| (?"^'^a/3 - Tar 13)+ r'^Sap 



(1.6) 



Here d is the dimensionality of space and _D is a con- 
stant characterizing the strength of velocity fluctuations. 
One assumes that the fluctuations are strong enough to 
ensure the large value of the Peclet number, that is. 



(1.7) 



The inequality ( |l.7| ) ensures the existence of the convec- 
tive interval of scales since it can be rewritten as L, 
where is the diffusive length 
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k/D. 



(1.8) 



The assumption of the Gaussian nature and zero cor- 
relation time for the flelds v and (f) allows one to derive a 
closed partial differential equation for the nth-order cor- 
relation function F„ of 6 ||l^,p8|,E7j . For the simultaneous 
pair correlation function F2(ri2) — {9{t,ri)0{t,r2)) one 
can solve the equation and find the explicit expression 
for F2. In the convective interval 



^2(r) = 2[F2(0)-F2(r)]~^: 



(1.9) 
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Comparing Eq. (1.9) with Eq. (1.3), one concludes that 
the exponent 7 introduced by Eq. (1^) directly deter- 



mines the scaling of the second-order structure function 
C2 = 7- 

However, for n > 2 the equations for Fn are too compli- 
cated to be integrated exactly. In |p^|-|l9|] the equations 
were analyzed in the limits 2 — 7 <gC 1 and ^7 3> 1, where 
the statistics of the passive scalar is close to Gaussian. 
The analysis led to an anomalous scaling that can be 
expressed in terms of the exponents Cn of the structure 
functions (1.2) and (1.3) 



Cn 



77.7 
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2(rf-|-2) 



n{n — 2) . 



(1.10) 



This expression covers both limit cases 2 — 7 <C 1 and 
drf ^ 1. The first term on the right-hand side of Eq. 
( 1.10 ) represents the normal scaling whereas the second 
one is just the anomal ous s caling exponent. The calcula- 
tions leading to Eq. ( 1.10 ) are correct if the anomalous 
contribution is much smaller than the normal one, which 
implies the inequality 



n < 



2^ 



(1.11) 



Below we will develop a different approach to the prob- 
lem. It will allow us to find the exponents C„ [Eq. (1.3)] 
of the structure correlation functions (1^) for any order 
n ^ 1 under the same additional condition ^7 3> 1 as in 
0,0- 



B. Path integral 

Generally, the statistics of classical fields in the pres- 
ence of random forces can be examined with the help of 
the field technique formulated in 1 29-31 1. In the frame- 
work of the technique, correlation functions are calcu- 
lated as path integrals with the weight exp(iX), where 2 
is the effective action related to dynamical equations for 
the fields. For the passive scalar in the Kraichnan model 
the effective action is 



ile^i J dt dr [pdt9 + pv\/e + KVpVS] 

-i J dtdridr2x{\ri - r2\)p{t,ri)p{t,r2) , (1.12) 

where p is an auxiliary field conju gated to 6. The first 
term in the effective action ( |1.12 ) is directly related to 
the left-hand side of Eq. (LI). The quadratic in the p 
term in Eq. ( 1.12| ) appears as a result of averaging over 
the statistics of the pumping (f>. 

Simultaneous correlation functions of 9 can be repre- 
sented as functional derivatives of the generating func- 
tional 



Z{X) = ( exp 



i / drX{r)e{t = 0,r) 



(1.13) 



where angular brackets designate averaging ove r the 
statistics of </) and v. With the help of the action ( 1.12| ) 
the generating functional can be rewritten as the path 
integral 



Z(A) = J VeVpVv cxp -J'iv) 
+ile+i J drX{r)e{t = 0,r) 



(1.14) 



Here y^{v) determines the statistics of the velocity field. 
Since we assume the Gaussian nature of the statistics, 
T{v) is a functional of second order over v with the 
kernel determined by the pair correlation function ( |1.5| ). 
Knowing Z{X) one can restore the probability distribu- 
tion function (PDF) of 9. It is convenient to treat the 
PDF of a particular object 



i9 



drf3{r)e{t = 0,r), 



(1.15) 



with a given functio n /3 (r). For example, the set of the 
structure functions (1.2) can be assembled into the PDF 
of the scalar differen ce in two points 9{r/2) — d{—r/2) 
which is the object ( 1.15 ) with /3(ri) — S{ri — t/2) — 
S{ri + r/2). The PDF of i9 is written as 



/ ^exp(-*yi?)Z[y/3(r) 

J-00 

Moments of ^ are then expressed as 

di? iz?rT'(i?) . 



(1.16) 



(1.17) 



We will be interested in the high-order correlation 
functions of 1!} or, in other words, we consider the limit 
n S> 1. This is eq uivalent to examinin g th e large tail 
of the PDF ( 1.16 ). One could expect [g^ that the tail 
can be calculated in the saddle-point approximation since 
there is a large parameter in the corresponding path 
integral. Unfort unat ely, direct application of th e method 
to the integral ( |l.l6 ) or to the moments ( 1.17 ) does not 
lead to success. 

To recognize the reason, let us consider the transfor- 
mation of the variables |3^] (see also ||3^ ) 

V ^ Xv , p~^Xp, t ^ X^^t , y Xy , k~*Xk. 

One can check that under this transformation all the 
terms in the square brackets in the right-hand side of 
Eq. (1.14) acquire the factor X, which means that in the 
saddle-point approximation In [2(?//3)] = yf{y/K) with 
some unknown function /. On the other hand, we ex- 
pect that correlation functions of the scalar itself (but 
not of its gradient, for example) do not depend on the 
diffusivity and the results of the works |15|-[l^,|2l|,p^,p8[ 
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confirm the expectation. Then, at small k the function 
/ remains a K-independent constant and we obtain 

lnZ(y/3)(x (1.18) 

Unfortunately, Eq. ( [l.l^ ) does not help to restore V{'d) 
since after substituting it into Eq. (1.16) we realize that 
the characteristic value of y in the integral can be es- 
timated as y ~ . Therefore, at large d the main 
contributio n to t he integral is determined by the region 
where Eq. ( 1.18 ) does not work. 

We conclude that the naive instantonic approach to the 
problem fails. The reason is that for the instanton the ve- 
locity field is fixed (does not fluctuate) in time and space. 
Obviously, a saddle-point solution is anisotropic because 
of the incompressibility condition drvv = 0. Fluctua- 
tions related to smooth variations of the anisotropy axis 
in time and space are strong and destroy the saddle-point 
approximation for the tail of the PDF Vid) or for the 
high moments of i?. Thus we should transform the prob- 
lem to more adequate variables, where fluctuations of the 
velocity are partly taken into account. This is the only 
chance to construct an instanton with weak fluctuations 
on its background. This is the goal of the next section. 



II. LAGRANGE FORMULATION 

As we mentioned above, the diffusivity n does not enter 
the result for the structure functions. Therefore, we will 
assume k — Q 'm all the following calculations. However, 
one should be careful since in this case it is impossible to 
deal with point objects. To provide a regularization, we 
should assume that the characteristic scales of the func- 



tion [3 in Eq. (1.15) are larger than rd- In addition, the 
scales are to be much smaller than L since we are going to 
examine correlation functions in the convective interval. 



In the diffusionless case the left-hand side of Eq. (1.1) 
describes the field moving together with the fiuid. Then 
it is natural to pass into the Lagrangian frame where the 
process is trivial. For that purpose we introduce La- 
grangian trajectories Q{t) that obey the equation 



dtg^v{t,Q). 



(2.1) 



We will label the trajectories by the positions of fluid par- 
ticles at t = 0: Q{t = 0) = r. Equation ( p_.l| ) (where k is 
omitted) can easily be solved in terms of the Lagrangian 
trajectories 



0(0, r) 



(2.2) 



Since we are interested in the field 6 at t = 0, due to 
causality the integration is performed over negative time. 
Therefore, Eq. (2.1) should be solved backwards in time. 

A simultaneous nth-order correlation function of 9 can 
be written as the product of n integrals ( ^.2D , averaged 
over the statistics of v and (p. In this representation, av- 
eraging over the pumping is very simple. For example, 
the two-point correlation function is 



dt (x(i?i2)). , (2.3) 

— OO 

i?i2(i) = i?(i,ri,r2) = \Q{t,r^) - Q{t,r2)\ . (2.4) 

The angular brackets ()„ in Eq. (^) denote averag- 
ing over the statistics of v only, since the statistics of 
is already accounted for there. Similar formulas can be 
written for correlation functions of higher orders. Once 
this is done , one can assemble them into the generating 
functional (1.13) 



Z(A) 



exp|-iydMridr2x(-Ri2)AiA2|^ . (2.5) 



where Ai,2 — A(ri_2). Calculating the moments of the 
object ( 1.15| ) in accordance with Eqs. ( 1.16| ) and (1.17) 
we get 



(l^l"> = 



AyM 

2lT 



(exp(-.FA-*2/^? + nln|79|))^ , (2.6) 



^x = Y I dtdridr2x(i?i2)/3(ri)/3(r2). 



(2.7) 



At this point, we would like to stress the close connec- 
tion between the statistics of the passive scalar and that 
of L agra ngian trajectories which can be seen from 
Eq. (lH). 



A. Statistics of Lagrangian separations 



Equations (2.5,^^6) show that the correlation functions 
we are interested in us are expressed via the a vera ge of 
exp(— JF\) over the velocity. Note that !F\ [Eq. (2/7)] de- 
pends only on the absolute values i?i2(t) of Lagrangian 
differences (^|J) . Therefore, instead of averaging over the 
statistics of v one could find the answer by averaging over 
the statistics of the Lagrangian separations R12 ■ Due to 
zero correlation time of the velocity field, the statistical 
properties of the field R12 appear to be relatively simple. 

To establish the statistics of R12 we start from the re- 
lation 

J-^dtRj2 - C12 = Rj2^Rl24vio. - V2a) , (2.8) 



followi ng fr om Eqs. (2.1) and (p.4|). As shown in Ap- 
pendix A 1. the average value of C12 is nonzero: 



{Ci2) = -D. 



(2.9) 



Next, exploiting the expression (1.5) for the velocity cor- 
relation function, one can find the irreducible pair corre- 
lation function 

2D 

((Cl2(<l)C34(t2))) = -rQl2.34S{h - t2) • (2.10) 

a 

The explicit expression for the function Q is rather cum- 
bersome 
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n2,34 



d+ 1 - 7 
4(d-l) 



^12 -^^34 



8(d- 1) 



34 



^7 (^12 



^13 ^ -^23) 



[ -R23 + ^1 



R 



p2 p2 



13 



p2 p2 



''23 



1 

■"14 



12 



14 



^24) 



^34 - -R13) + (-^12 + ^24 - R14) {R. 

"■24 



■34 



p2 p2 



23) 



R 



34 



p2 p2 

^23 -"-24 



(2.11) 



It can b e fo und from the definition of C12 [Eq. (2.8)] 
formula (1.6) and relations such as 



-R12-R1 



kRi2 



-R23) : 



■R12-R34 ^2^^ 



^23 -"-13 



-R24) 



In the spirit of the conventional procedure [p9|-|31[ , one 
can assert that any average over the statistics of R12 can 
be found as the path integral over R12 and over an aux- 
ihary field mi2 = m{t, ri,r2) with the weight 



exp 



[ j AtAriAr2 (mi27 ^9ti?72 " '^i2Ci2) 



where angular brackets mean averaging over the statistics 
of the velocity. Since C12 is 5 correlated in time, the av- 
erage can be expressed in terms of Eqs. (2.9) and (2.1C) 
only. The result is exp {Hr) , where 



Hr = i J dt J dri dr2 mi2(7^^9ti?72 + -D) 



dri dr2 dra dr4 (3i2,34?™i2m34 . (2.12) 



Now, instead of Eq. (2^) we can write 

(1^1") =y J Pi?Pme''^«-^--'^^+"i"l''l. (2.13) 



The integration in Eq. (2.13) is performed over func- 
tions of t, ri, and r2 with some boundary conditions 
imposed on them. The condition for the field R12 follows 
from ^(0) = r and reads 



(2.14) 



The boundary condition for the field mi2 should be 
mi2(— 00) = 0, since we deal with free integration over 
R12 in the remote past. Note that due to the definition 



(2.4) the triangle inequalities 

i?12 + -R23 > R 



13 



(2.15) 



should be satisfied for any three points. Actually, the in- 
equalities are constraints that sho uld b e imposed on the 
field i?i2 when integrating in Eq. (2.13). 



B. General instantonic equations 



In the preceding subsection we derive a formula ( 2.13| ) 
for di^l"). Its calculation is equivalent to solving some 
nonlinear field theory. It looks very infeasible to perform 
this task. We are going to calculate the integral (2.12) 
in the saddle-point approximation regarding the number 
71 large enough. To be consistent, when doing the proce- 
dure one should remember about the constraints (2.15). 
Unfortunately, it is very hard to take them into account 
explicitly. We will ignore the constraints, which is cor- 
rect under the following conditions. First, the inequal- 
ities (2.15) should be valid in the instantonic solution. 
Second, fluctuations on the background of the instanton 
should be weak (this is also the applicability condition 
of the in stan toni c for malism itself). We argue in Ap- 
pendixes B 1 and B 2 that those conditions are satisfied 
if 



(i7 > 1 



(2.16) 



Note also that for the condition ( p. 16 ) fluctuations of 
a Lagrangian separation near its average value are weak 
(see Appendix A 2). The inequality (2.16) will be implied 
below. 



Thus we obtain from the integral (2.13) in the saddle- 
point approximation 



(2.17) 



Here solutions of the instantonic equations should be sub- 
stituted which are extrema conditions for the argument 
of the exponent on the right-hand side of Eq. (2.13). 
Variation over mi2 and R12 gives the instantonic equa- 
tions 



D f 

i{^^^dtR'l2 + ^"^ / ^''S ^''4 '312,34^-34 : 



(2.18) 



iR\2 ^dtmi2 + 



D 



dr3dr4 < 2 — mi2m34 



dRi2 



. C?Ql3,24 

-4^TT^ 7ni3m24 



y 



X'(i?i2)/3(ri)/?(r2). (2.19) 



dRi2 ■■ -^■■ -] 2 
The extremum conditions over y and -d read 



d^iy j didridr2x(i?i2)/3(ri)/3(r2), (2.20) 
iy^nl-d. (2.21) 
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Note that only Eqs. ( ^.18 ) and are true dynamical 

equations, carrying the i nform ation abou t the dynamics 
of the flow, whereas Eqs. (2.20) and (2.21) are constraints 
impos ed on the in stant onic solution. One needs to add to 
Eqs. (2.18) and (2.19) some boundary conditi ons. The 
value of the field R12 is fixed at t = by Eq. ( 2.14 ). As 
for the field mi2, we already noted that it should tend 
to zero when t —00. It can be understood as the 
extremum condition that appears after variation of the 
effective action over the boundary value of R12 in the 
remote past. 

One can easily establish the asymptotic behavior of 
R12 at \t\ 00. There the field R12 grows and loses its 
dependence on ri 2- The field mi2 tends to its vacuum 
zero value at \t\ — >_oo. Therefore, at large \t\ the term 



with mi2 in Eq. (2.18) can be omitted and we find 



i?T « -fD\t\ 



(2.22) 



The expression ( 2.22 ) is nothing but the Richardson law 
for divergence of Lagrangian trajectories Let us 

stress that now it holds on the classical (mean-field) level, 
without taking into account fluctuations on the back- 
ground of the instanton. To clarify this point, notice 
that if the velocity field is a deterministic function of 
time and space (as it is for the naive instanton discussed 
above) then the Richardson law can not be valid for all 
the Lagrangian trajectories. In our instanton we get rid 
of the velocity field that resulted in the emergence of 
t he R ichardson law. Note that the tri angle inequalities 
( ^.15 ) are obviously satisfied both for ( 2.14 ) and for the 

asymptotic behavior (2.22). 

The expression for the action appearing in Eq. ( 2.13| ) 



iX = iXR-Tx = i J dtdridr2-f~'^mi2dtRj2-E , (2.23) 
^ " y/ d^id^2X(^i2)/3(r-i)/3(r2) - iD J dridr2mi2 



/ dri dr2 drs dr4 (5i2,34mi2m34 . 



(2.24) 



We see from Eq. ( 2.23| ) that the quantity E plays the 
r ole o f the H amil tonian function of the system, while Eqs. 
(2.18) and (2.19) arc canonical equations corresponding 
to the Hamiltonian function. Since E does not explicitly 
depend on time t, its value (which can be called energy) 
is conserved. Actually the energy is zero on the instan- 
tonic solution since at t — )■ — 00 we have mi2 and 
R12 00. Note that since the Hamiltonian ( 2.24 ) ex- 
plicitly depends on the coordinates via /3, there is no 
momentum conservation law. 

Before proceeding to the solution of the instanton 
equations, let us make a remark concerning fluctuations 
on the background of the instanton. In the linear ap- 
proximation over the fluctuations we obtain an estimate 
for the typical fluctuation of 



Note that the fluctuations of R tend to zero when t ^ 
since R12 i s fixe d at t = 0. Comparing the estimate ( 2.25 ) 
with Eq. ( ^.22D we obtain 



{5R^f/R^^ - d- 



(2.26) 



We conclude that the fluctuations on the background 
of our instanton are weak provided d 3> 1. The above 
evaluations are rou gh a nd need a more accurate anal- 
ysis (see Appendix B2). Neve rtheless, they show that 
the Richardson behavior ( 2.22 ) inherent for our instan- 
ton suppresses fluctuations on its background. 

The system ( ^.18 ) and (2.19) consists of two nonlinear 
integro-differential equations with boundary conditions 
imposed on the opposite sides of the time interval, that 
is, at i = for R12 and at i = —00 for mi2. Therefore, in 
the general case it is very difficult to solve the instanton 
equations. Nevertheless, one can hope that for some par- 
ticular objects the system of equations can be reduced 
to a simpler form allowing the complete solution. This 
hope comes true for the structure functions. 



III. INSTANTON FOR STRUCTURE 
FUNCTIONS 



Using the general scheme developed in Sec. || we 
will examine the expressions for the structure functions 
(1.2) at large n. In other words, we will be interested 
in the statistics of the passive scalar difference taken at 
the points separated by the distance r. Since the dif- 
fusivity is neglected, we cannot examine the difference 
d{r/2) — d{—r/2) itself. Nevertheless, we can treat the 
statistics of the differences averaged o ver se parations near 
r. So we should consider the object ( 1.15 ) with 



(3.1) 



Here 6 Air) is the function with the width A^^ ^ sat- 
isfying the condition J dr 6\{r) = 1, which can be called 
a smeared S function. Then we can write 

5„«(m~exp(zX-n + nln|^|)li„,, , (3.2) 



where we used Eq. ( 2.17 ) and substituted Eq. ( 2.21 ) 



A. Reduction 



{SR^y r^-fDR^\t\d- 



{2.25) 



Now we turn to the instantonic equations ( 2.18| ) and 
( ^.19 ). Let us observe, that since the source on the right- 
hand side of Eq. ( 2.1S| ) is proportional to /3(ri)/3(r2), the 
field mi2 can be approximated as 
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mi2 



im^ ( / r\ / r 



(3.3) 



where m± are functions of time only. Writing it, we im- 
plicitly assumed that the field R12 is smooth near the 



points ±r/2. Then the relations (2.20) and (2.21) give 



2n J dt {x{R+)-x{R~)} , 



(3.4) 



where we introduced 

R+it)=R{t,r/2,r/2), R^{t) ^ R{t,r /2,-r/2) . (3.5) 



Sub stituting the expression (|3.3| ) into Eqs. ( 2.18 ) and 
(2.19) we obtain a closed system of ordinary differential 
equations for m± and R± . It is convenient to proceed in 
term s of t he effective action. Substituting Eq. ( |3.3| ) into 
Eq. (|2.23D we get 



.0 

il= dt [7-1 {m^dtRl - m+dtRX) ~ ^] 

J — OC 

E = y^x {R+) - X {R-)} + D{m+ - m_) 

£'(2-7) ,9 2 1 

^ Ad{d- 1) y^-^^ + 2m-m+ip2 + m^ip^i] . 



(3.6) 



(3.7) 



Here we introduced the designations 



VI = — 7, R~ 

2-7 



R 



2-7 



i?1 



[Rl - Rl] 



27-4 



R 



(P2 = R 



R 



4-7 ^ (2fl^-fl^)^^ 



R 



2-7 



2i?2^ - R\ 



(3.8) 



2-7 



^3 = -i?^ 



Since the effective action ( p.6[ ) depends only on the func- 
tions TO± (t) and i?± {£) , one can obtain the system of 
ordinary differential equations for the functions as ex- 
tremum conditions of the action. The boundary condi- 
tions for the e quations are i?+ = and i?_ = r at t = 
[see Eq. (2.14)] and m± ^ at < — > —00. Resolution of 
the system allows one to find m± and R± as functions of 
time. Once they are known, it is possible to resto re the 
function R12 in the whole space from Eq. ( 2.18 ). The 
problem is discussed in Appendix B 1 . There we argue 



that the function i?i2 is really smooth in space which is 
a justification of the pr ocedu re described. 

Since we accept Eq. ( ^.16 ), d 3> 1. Using the inequal- 



ity one can keep in the functions (3.8) only the terms of 



the main order over d. This means that one can neglect 



in Eq. (3.S) the second contribution to (pi in compari- 
son to the first one and also </32, V3 in comparison to Lp\. 
Potentially this procedure is dangerous. We will show 
that due to the smallness of r/L, the intervals where 
i?„ — i?+ ^ i?_ play an important role. Then, we see 
that it is the difference of R± that enters the first term 
in iy9i, while the others do not contain this smallness. 
Therefore, we observe cancellations that could lead to a 
competition of d and Ljr (the latter parameter is con- 
sidered as the largest in the problem) . To check the pos- 
sibility, we performed calculations keeping all the terms 
The calculations are sketched in Appendix 



(3. 



in Eq. 

They show that in the final expressions only combi- 
nations of iy9i.2,3, containing the same cancellations are of 
importance. The legitimacy of the procedure is proved. 

Omitting '£23 in the expression (3.7) and then varying 
the action (3^) over m+, we get a trivial equation for 
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-^dtRl 



D. 



(3.9) 



Its solution, satisfying the boundary condition -R+(0) = 
is simply 



Rl=lD\t\ 



(3.10) 



To examine the behavior of _R_ it is convenient to pass 
to the new variables 



i?4 



fi = m^R\ . (3.11) 



As time t goes from to —00, the variable ^ runs from 
—00 to -l-oo and v runs from -l-oo to 0. The latter is 
clear from the asympto tic b ehavior _RZ w RI^ = ^D\t\ at 
t — > — 00. The relation (3.4) in terms of the new variables 
is 



-t-00 

2n^ J d^e^« [x{R+)~xiR-)] 



(3.12) 



Recall that the energy E entering the action (B.6) is 
an integral of motion whose value is equal to zero. Thus 
we can perform the standard procedure of excluding a 
degree of fr eed om in a canonical system. Equating the 
expression (3/7) to zero, we can express m+ in terms of 
/i, V, and ^. Substituting the result into Eq. (3.6), we 
get 



-\-oo 

-il= J d^{j-'fid^v-H), 



(3.13) 



ff=-M^; + ^0(«) + ^^[x(i?+)-x(i?-)]e^«, (3.14) 

^=(l+t;)2-4/7 r(l+^)2/7-l_ll r(l+„)2/7_ll (3,15) 



Here we kept main contributions over d only. I n En . 
( |3.14 ) we set = — |yp since as follows from Eqs. ( 2.21 ) 
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and (3.4) y is an im agin ary number. Extremum condi- 
tions for the action ( 3. IS ) read 



_i du 



dH 



dH 

dv 



1^^ (3-16) 



which are canonical equations f or th e variables and v 
in the time ^. Of course, Eqs. ( 3.16| ) could be obtained 
directly from the extrema conditions for the action (3.6). 

To conclude, we reformulated the problem as follows: 
find such a value of y that the solution o f Eqs . (3.16) with 
the given y, being substituted into Eq. ( 3.12| ), reproduces 
the correct value of "i? = —in/y. Below we discuss the 
first and the most diffic ult p art of the program that is 
solution of the system ( 3.16|) . Though it cannot be in- 
tegrated exactly, we can solve the system approximately 
by asymptotic matching, which is enough to determine 
the structure functions S„ . 



B. General structure of the instanton 

The evolution of i?_ in 'time' ^ can be divided into 
three stages. During the first stage, starting at ^ = -co, 
both i?+ and R- are much less than L and it is possible 
to substitute bot h and x(^-) by x(0)- Then the 

last term in Eq. (3.14) is equal to zero. Du ring t he second 
stage R± ~ L and the last term in Eq. (|3.14|) is of im- 
portance. During the final stage, where ^ L, 
one can again neglect the last term in Eq. ( ^.14 ). Note 
that only the seco nd sta ge contributes to d'^, which can 
be se en from Eq. ( 3.12 ). Since the Hamiltonian H [Eq. 
(3.14)] does not explicitly depend on 'time' ^ during the 
first and the third stages, its value is conserved there. 
Actually, the value of H is equal to zero during the third 
stage since fi ^ and R+ « R- :s> L a.t ^ +oo. On 
the other hand, during the first stage the value Hi of the 
Hamiltonian function H is nonzero. Therefore, during 
the second stage the value of H diminishes and should 
finally reach zero when the trivial third stage starts. The 
value of Hi as a function of n has to be established from 
the matching of the stages. 



Now we are going to solve Eq. ( |3.16| ) for the first stage. 
Resolving the equation H = Hi in terms of /x we get 



Giv)=±\ + 



iHi^ 



Then we find from Eq. ( |3.16| ) 



(3.17) 
(3.18) 

(3.19) 



At ^ ^ — oo (that is at small \t\) the function v should de- 
crease with increasing ^ since i?_ « r and _R-|_ incre ases. 
To ensure the negative value of dv/dS, in Eq. (3.19) one 



shoul d take the positive sign of the square root in Eq. 
(|3.18 ) , which leads to a positive value of G. The sign of 
G can be changed if during the evolution G turns into 
zero which corresponds to the presence of a reverse point 
in the dependence v on ^. 

Equation (3.19) enables one to find t; as a function of 
^. Let us integrate the equation over ^ from — oo to some 
value. Then we get 



da: 



G{x) 



In 



vRl 



(3.20) 



To avoid difficulties related to infinite values of ^ and v at 
the initial point, we subtracted from G^^ its asymptote 
G~^{x) « 1/a; a t larg e x. This enforces the convergence 
of the integral (3jC) at large x. The constant of in- 
tegration in Eq. (|3.20 ) was established from the limit 
w — > oo: since the integral on the left-hand side of Eq. 



( 3.20 ) tend s to zero as v increases, the right-hand side of 
Eq. (3.20) should also tend to zero. This r equir ement 
is ensured by the r-dependent factor in Eq. (3.20) since 



iVj^ ~ r'^ / V at V oo, as fol lows from the boundary con- 



dition i?_(0) = r and Eq. ( |3.1l| ). The left-hand side of 
Eq. (3.20) should be viewed as a contour integral, which 
determines its value in the case of the non-monotonic 
behavior of f as a function of ^. 

Equation ( 3.2C ) allows us to establish a relation for the 
parameters characterizing the first stage. Let us consider 
the integral over the wh ole fi rst stage. Then we should 
substitute v = v^, in Eq. (3.2C), where is the value of v 
at the end of the first stage. The initial substage (where 
w ^ 1) gives a constant o f ord er unity in the integral on 
the left-hand side of Eq. ( 3.20| ) since G{x) « x there. We 
neglect the contribution substituting u ~ 1 as the lower 
limit in the integral. Then the integral J dx/x produces 
just Inu, which is canceled by the corresponding term on 
the right-hand side. Next, at the boundary between the 
first and the second stages i?_ ~ L since the pumping 
enters the game there. Therefore, with the logarithmic 
accuracy one can write 



dx 



7 In 



L 



(3.21) 



We see that there is a large parameter L/r in the ar- 
gument of the logarithm on the right-hand side of Eq. 
( |3.21 ) . An analysis shows that due to this large parame- 
t er th ere are only two possibilities to satisfy the relation 
( |3.2lD . Both of them are related to zeros of the function 



G because only near the points where G is small can the 
integral reach a large value. The first possibility is real- 
ized when G is zero only at v — 0. In this case <C 1 
and w is a monotonically decreasing function. The sec- 
ond possibility is that G has a zero at some point v = Vr- 
That is just the reverse point where the derivative dv/d(, 
changes its sign; see Eq. (3.19). Then, the integral on the 
left-hand side of Eq. ( 3.21[ ) is determined by the vicinity 
of the point since G is small there. 
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A choice between the possibihties depend on the value 
of Hi. If Hi > He, where 



Hr 



8(2-7) 



then G cannot be zero (except for the point v 
the integral on the left-hand side of Eq 



(3.22) 



0). Then 
(3.2l|) reaches 



3^) 



its large value at u ^ 1. Substituting into Eq. ( 
asymptotic expression 

2(2-7) 2 



the 



(3.23) 



valid at ?; <C 1, we can c alculate the integral on the left- 
hand side of Eq. (3.21) with the logarithmic accuracy 
and find 



Inu* = "f^/T 



Hi I He In f^- 



(3.24) 



We see that due to r <C i, indeed ^ 1. 

In the opposite case Hi < He the situation is more 
complicated. From the asymptotic expression 



27 



(3.25) 



valid at w <C 1, we see that G is zero at v — Vr, where 



16(2-7) 
d7(4 - 7) 



(He -Hi). 



(3.26) 



It is just the reverse point where the derivative dv/dS, 
changes its sign. Therefore, the sign of G is positive dur- 
ing the initial part of the first period and negative during 
t he fi nal one. Thus we should take the upper sign in Eq. 
(3.18) for the first part and the lower sign for the second 
part. The main contribution to the left-hand side of Eq. 
(3.21) is determined by the region near the re verse point 
V ~ Vr ^ Vr where we can use the expression (3.25). The 
explicit integration gives 



2(2 - 7) VHc - Hi 



: 7 In — . 

r 



(3.27) 



Since the logarithm is large, Hi is close to He and hence 
<C 1, as we im plicit ly assumed in the expression (3.25). 
Note that Eq. ( |3.27 ) does not fix the value of w*, as it 
was for Hi > He- 

Now, we shou ld extra ct add itional relations that, along 
with Eq. ( 3.24 ) or Eq. ( 3.27 ) will fix the instantonic so- 
lution and determine the final answer for the structure 
functions. It can be done by establishing the evolution 
during the second stage and by its subsequent match- 
ing with the first stage. Unfortunately, the procedure is 
rather lengthy and is individual for each particular case 
We present the calculations in Appendix 



C. Expressions for structure functions 

Based on the reasoning given in the preceding subsec- 
tion and on the calculations described in Appendix 
one can establish exp ressions for the structure functions 
from the relation (3.2). Here we enumerate basic results, 
referring the reader interested in technical details to Ap- 
pendix O. 

The case Hi > He is realized ii n < Ue (see Appendix 



DlJ), where 



dj 



2(2-7)- 



(3.28) 



Calculating the action 1 and ■& (see Appendix D 1 ) and 
substituting the result into Eq. (3^), we obtain 



nP2Ci 



r 

J D " J \L 

(2-7)^2 



i/2 



2d 



(3.29) 
(3.30) 



The quantity Ci in the expression ( |3.29| ) is a constant of 
order unity, whose value depends on the shape of x (that 
is on the details of the pumping) and is consequently 
nonuniversal. Note that the r-independent factor in Eq. 



29 ) is determined by the single-point root-mean-square 
value of the passive scalar 



(3.31) 



Comparing the expression ( 3.30 ) with Eq. ( 1.10| ), we see 
that they coincide under the conditions n 3> 1 and d 3> 1 
that were implied in our deriva tion. Surprisingly, the n 
dependence of C„ given by Eq. ( 1.10 ) is correct not only 
in the limit (1.11) (that is for n <ti Ue), but up to n = U e, 
which is the boundary value for Eqs. ( ^.29| ) and ( 3.30 ). 
A detailed consid erati on of the case Hi > He is pre- 
It shows that this possibility 
Then, the scaling exponents C„ 



sented in Appendix D 2 
is realized at n > ? 



appear to be n independent and equal to the value 



Cc 



d-f^ 



8(2-7)- 



(3.32) 



The n-dependent numerical factors in S'„ can be found in 
two limits: n — Uc <^ Uc and n ^ Ue- The former case is 
discussed below, while in the latter case one can obtain 
(see Appendix D2) 



nP2C2 
D ' 
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n/2 



(3.33) 



The quantity C2 in Eq. (|3|) is again a non-universal 
constant of order unity. The expression (3.33) corre- 
sponds to the factorized Gaussian PDF 
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exp 



(3.34) 



Let us stress tha t wh en calculating Sn ~ d"!?!") with the 
help of the PDF ( [3.34| ) , the characteristic d is of the order 
of the s ingle -point root-mean-square value of the passive 
scalar ( ^.31 ) and the relatively small value of the result 
( |3.33 ) compared to a single-point val ue is e nsured only by 
the small r-dependent factor in Eq. ( |3.34 ). In Appendix 
we establish the inequality 



n L 
in — < 7 In — 

a r 



(3.35) 



which restricts the region where expression (3.33) is cor- 
rect. For larger n the character of the PDF essentially 
changes and it te nds to a single-point PDF that is sim- 
ilar to Eq. ( 3.34 ) but does not contain the r-dependent 
factor. 

Note that the cases 7 <C 1 and 2 — 7 <C 1 need a special 
analy sis which is performed in Appendix p 3| . The answer 
( |3.33 ) should be slightly corrected in the case 7^1 and 
keeps its form at 2 — 7 <C 1. 

We can treat the structure function Sn as a continu- 
ous function of n. Then the vicinity of the critical value 
n = ric requires a s epar ate consideration which is pre- 
sented in Appendix D4. The main peculiarity that ap- 
pears in the expressions for the structure functions is a 
critical dependence on n. The expression for the struc- 
ture functions can be written as 



in 



iric 



D 



a/2 



(3.36) 



which implies the condition |n — nd <C tIc- The factors 
C± are non-universal constants of order unity which are 
different for the c ases n < and n > Uc- The ex ponen ts 
(^n in expression ( ^.36| ) are determined by Eq. ( ^.30 ) if 
n < Uc and Cn = Cc [Eq. ( ^.32| )] \i n > ric- In the 
consideration made above we suggested that r/L is the 
smallest parameter of our theory. However, if n — s- ric, 
then jric — n\ starts to compete with r/L and at small 
enough ric — n the consideration presented in Appendix 
D 4 is inapplicable. The criterion that determines the 



validity of Eq. (3.36) is established in Appendix D4 



L Uc 
7ln- > I 1 



(3.37) 



We see that the first factor in Eq. ( 3.36 ) possesses the 
critical behavior proportional |rt — nd"'' that is saturated 
i n the narrow vicinity near n = ric, where the condition 
( |3.37 ) is violated. To avoid a misunderstanding, let us 
stress that despite the critical behavior, S'„ remains a 
monotonically increasing function of n at a fixed L/r. 
This is obvious for n > ric, whereas for n < it ac- 
counts for the stronger dependence on n of the second 
(r-dependent) fact or in Eq. ( [3.36| ), which is guaranteed 
by the inequality (3.37). 



We presented the results of the analysis based on the 
saddle-point approximation. The account of fluctuations 
on the background of our instanton could, in principle, 
change the results. Particularly the value of (n could 
increase. Therefore, one should estimate the role of the 
fluctuations. The corresponding analysis is presented in 
Appendix B 2 . It shows that for the condition ( 2.16D fluc- 
tuation effects are weak and cannot essentially change the 
results obtained. 



CONCLUSION 

We have performed an investigation of the structure 
functions in the Kraichnan model in the framework of 
the instantonic formalism. Though our approach is cor- 
rect only for large dimensionalities of space, we observe 
a nontrivial picture, some peculiarities of which could be 
realized in a wider context. Below we discuss the results 
obtained. 

We have established the n dependence of the scaling 
exponents, which are determined by the expr ession ( 3.30| ) 
for n < ric and remain the constant (3.32) for n > ric 
where ric is defined by Eq. ( [?.2^). Our results contra- 
dict the schemes proposed in p6| , |27| | . The value ( ^.32 ) is 
different from and smaller than the constant obtained in 
psf , which can be considered really as an estimate from 
above. For n ric our expression coincides with the 
answer obtained perturbatively (l^Jl^ at large d. Sur- 
prisingly, the quadratic dependence of C„ on n is kept up 
to n = ric- Such an n dependence of ^„ is well known 
from the so-called log-normal distribution proposed by 
Kolmogorov pO| . 

The expressions (3.2S) and ( ^.33| ) reveal the combina- 
toric prefactors in Sn that are characteristic rather of 
a Gaussian distribution. A natural explanation can be 
found in terms of zero mode ideology |p"5 Hl9| , ^ . We 
know that for n > 2 the main contribution to the struc- 
ture function Sn in the convective interval is related to 
zero modes of the equation for the nth-order correla- 
tion function of the passive scalar. The exponents of 
the modes are determined by the equation (and could 
be very sensitive to the value of n), whereas numerical 
coefficients before the modes (determining their contri- 
bution to Sn) have to be extracted from matching on the 
pumping scale where the statistics of the passive scalar 
is nearly Gaussia n. Th is exp lains the combinatoric pref- 
actors in Eqs. (|3.29 ) and ( 3.33 ). Probably the most 
striking feature of our results is the unusual behavior of 
Sn (treated as continuous functions of n) near n = ric. 



which is determined by the expression (3.36). 

Now we briefly discus s the interp retat ion of our results. 
The log-normal answer ( |3^ ) and ( |3!30| ) can be obtained 
if we accept that for large fluctuations, giving the main 
contribution into the structure function Sn , the pumping 
is inessential and the fluctuation is smooth on the scale 



Then one obtains from Eq. (1.1) the equation for the 
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passive scalar difference taken at the separation r, 

dt\n{M) = -vr/r'^ , (3.38) 

where we substituted V6' by {A6)/r. We immediately get 
from this equation a log-normal statistics for A9 that is 
a consequence of the central limiting theorem. The sat- 
uration at n > ric can be explained by the presence of 
quasidiscontinuous structures in the field 6 making the 
main contribution to the high-order correlation functions 
of 6. Note also a similar non-analytical behavior of Cn for 
Burgers's turbulence 0, which is explained by the pres- 
ence of shocks in the velocity field. Although formally 
our scheme is applicable only in the limit ^7 3> 1, one 
can hope that the main features of our results persist 



for arbitrary values of the parameters. This hope is sup- 
ported by [ [l3[ , where a saturation of was observed in 
numerical simulations of the Kraichnan model at d = 'i. 
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APPENDIX A: SINGLE LAGRANGIAN 
SEPARATION 



where R = R{t). 
Let us write now 



In this appendix we treat the statistics of a single La- 
grangian separation defined by Eq. (2.4). The consider- 
ation will allow us to est ablish the relation (2.9) and also 
to clarify the condition (2.16). 



t-At 



1. Richardson law 

A single Lagrangian difference R between two La- 
grangian trajectories g and ^ -f- i? is governed by the 
equation 

dtRa = Wa ^ Vaig + R) - Vaig) , (Al) 
with the correlation function 

{wc,{ti)wf3{t2)) = 2/C«/3(i2)(5(ti - t2) , (A2) 
IC^p{R)^^R-^ !^-^^^{R^Sai3 - RaRp)+R^S^p 



follo wing from Eq. (|L6|). First of all we get from Eqs. 
(^.10y2.1l|) 



(C(il)C(i2)) = ^i?^<5(tl-i2) 



(A3) 



where in accordance with (2^) C = 7 ^dtK^ . Then 

t-At 

R{t- At)- R{t)^ R^-''{t) I dTC(T), 



(A4) 



where we believe At > to be a small time interval 
(remind that we treat evolution backwards in time). Av- 
eraging over the velocity sta tistics o n the interval from t 
tot - At we get from Eqs. (|A3p^ 



9n 

[R{t-At)-RY) = ^^'"'^ 



(A5) 



Rc,[t - At) ^ R^[t) + j dTu;„[r,i?„(T)], (A6) 



which is the direct consequence of Eq. (Al). Then we 
find from Eq. (|A6D in the approximation needed for us 



f-Af 



R^{t-At)-R^{t)K2Ra j dTWa{T,R) 



t-At 



dr dr' Wa (t, R) Wq (t' , R) 



t-At 



+i?Q^^ J dTWa{T,R) j dT'wpir' ,R) , 



where again R = R{t) and we used the incompressibility 
condition dw/dR = 0. Averaging the expression over 
the velocity statistics on the interval from t to t — At we 
get 

(i?2(< - At) - w 2^{d + 2 - -i)R^-^At , (A7) 



where we used the expressions (A2) and taken into ac- 
count dlCap{R)/dRa = 0. Then we obtain from Eqs. 



(R(t-At)-R)^^(d+l--f)R^-^At. (A8) 
d 



(A9) 



Next, we get from the definition (|2^ 

{R^{t-At)-R'') =-7(0 At. 
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Expanding here the difference up to the second order 
over R{t — At) — R and substituting then ( A5,|A^) we 
find finally 



(0 



-D. 



(AlO) 



Note that the average is negative which is a consequence 
of considering an evolution backwards in time. 

The average value (C) is obviously the same for al l the 
Lagrangian separatio ns. Th erefore, we arrive at Eq. (2.£ ) 
leading then to Eq. ( 2.22 ), which is a manifestation of 
the Richardson law. 



2. Simultaneous PDF 

The Kraichnan model admits a closed description of 
the simultaneous probability distribution function V{R) 
for a ny s ingle Lagrangian difference R. The point is that 
Eq. (Al) in this case can be considered as a stochastic 
process with white noise on the right-hand side. It is 
well known how to obtain the e quat ion for 'P{R) in the 
situation. Using the expression (A2) we get 



dV/d\t\^lC^p{R)V^VpV . 



(All) 



The same problem for 7 = 2/3 was considered by Kraich- 
nan who obtained the asymptotic behavior of V at 
large times t. Here we present a straightforward general- 
ization of Kraichnan's scheme for arbitrary 7. 

Suppose that at t = the PDF is "P = Rl~'^d{R- Rq), 
the initial condition corresponds to a fixed value of \R\ 
at i = 0. The above expression implies the normalization 
condition J dR R'^^^'P{t, R) — 1. Solutions correspond- 
ing to other initial conditions can be exp resse d via this 
fundamental solution since the equation (A12) is linear. 
Due to iso tropy , V will be a function of R only. Then the 
equation (All) is rewritten as 



Dd\t\ Rd-^dR\ OR 



(A12) 



Of c ours e (A12) can also be obtained directly from Eqs. 

([a^,|aio|). 

Performing Laplace transform of Eq. ( A12| ), one ob- 
tains 



V{t,R) 



^^e^^n^X\t\]S{\,R), (A13) 



A—ioo 



where all the singularities of 5(A) have to b e to t he left of 
the integration contour. The function S in (A13) satisfies 
the equation 



1 



d 



Rd-^ dR 



dR / Rt 



j5{R~Ra). 



The equation can be solved separately in the regions 
R < Rq and R > Rq where we deal with the homo- 
geneous equation and then the matching conditions at 



R = Rq give us coefficients. Assuming suitable bound- 
ary conditions (5 is finite at i? — > and at i? — > 00) we 

get 



5(A,i?) = -(i?i?o)^'/'+^/' 
7 



(A14) 



U {^R''^') kJ^R^,^'] if R<Ro, 
(2f J. if i?>i?„. 

Here, v — d/j — 1 and lu is the modified Bessel function, 
is the McDonald function. 
If we are interested in the asymptotic behavior of 
V{t,R) at large times \t\ > dR'^/{D-/'^) then R near the 
maximum of V satisfy i? S> i?o- In addition, in the inte- 

cx is small. 



cha 



gral ( |A13| ) the characteristic value A, 
Thus we can take only the second term in (A14) (corre- 
sponding to R > Rq) and substitute the first term of the 
expan sion I,y{z) ~ r~^(l -I- z^)(z/2)'^. Then the integral 
(AI2) can be taken explicitly and we obtain 



rit,R)R'^-^dR 



1 



r(i + i.) 

dR'' 



derexp(-0: 



Dj^\t\ ' 



(A15) 
(A16) 



which gives the asymptotic self-similar behavior of PDF. 
We see that \nV oc —R''. It is interesting that the asymp- 
totic PDF is Gaussian at 7 = 2 for any d. As we move 
7 from two to zero, the PDF is getting more and more 
non-Gaussian reaching an extreme non-Gaussianity (log- 
normality) at 7 = 0. Note an obvious consequence of 
Eqs. (|A15|JA16D: 



(A17) 



for arbitrary /i. The asymptotic law ( A17 ) is a manifes- 
tation of the Richardson law. 

It is not very difficult to calculate first moments using 
the expression (A15): 



(R') = 



Df[t\ 
d 

Dl'\t\ 



2/7 



4/7 



r(i 



2/7) 



r(i 



r(i 



■ 4/7) 



r(i + i.) 



In the large d limit the expressions (A18,A19) give 



{R') = {Dj\t\r/\ 



{R'y 



exp 



id) 



(A18) 

(A19) 
e 

(A20) 
(A21) 



irrespective to the value of 7. We see that fluctuations of 
i?2 are small if 7^ ^ 1 and in the opposite limit 7^ <C 1 
the cumulant of R^ is much larger than (i?^)^. We con- 
clude that the inequality ( 2.16| ) is just the cond ition at 
which R'' only weakly fluctuates near the value ( ^.221 ). 
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APPENDIX B: 

Here we discuss the applicability c ondit ions of our 
scheme. First, the triangle inequality (2.15) should be 
satisfied for our instantonic solution. Second, fluctua- 
tions on the background of our instanton should be weak. 



1. Triangle inequality 



Here we discuss the triangle inequality ( 2.15| ) that was 
ignored in our saddle-point calculations. We ar gue that 
if the parameter c?7 is large, the inequality ( 2.15| ) is satis- 
fied for the instanton solution. In other words, we can say 



that calculating the path in tegra l (2.13), we can dismiss 
the restriction supplied by (^.15|) , because the contribu- 
tion from the regions where the inequality is violated, is 
small in this limit. 



Let us recall that the inequalities (2.1£) are obviously 
satisfied both for the initial condition (^.14 ) and the 
asymptotic behavior (2.22). Therefore, they could be 



violated only for times of the order of the instanton life- 
time. Our project will be as follows. First, we check that 



the inequality (2.15) holds for the instanton solution if 
we consider the instantonic equations in the main order 
over \/d. Then, we will show, that next order terms over 
\/d considered in Appendix |c| can lead to violation of the 
inequality if ^7 ^ 1. Generally, this procedure is compli- 
cated, and here we present only some calculations which 
serve as a basis for our conclusions. 

As a first step, we have to restore the field R12 in the 
whole space, that is for any two points ri a nd r2 . This 
can be done as follows. Le t us consider Eq. ( 2.18| ). Sub- 
stituting mi2 in the form (|3.3D we get 



l-^dtRl^ + D 



2d 



(Bl) 



rI-^+rI+^-rI+^-rI-'^j {RI-+RI+-RI+-R2 



We kept only the main term over 1/d m Q ( ^.ll| ) . In Eq. 
(Bl) we introduced auxiliary fields 



=i?(ri,r/2) , =i?(ri,-r/2) 
They satisfy the closed system of equations 



l-'dtRj+ + D = - 
l~^dtRj_+D^- 



Rl-J + R^^ - Rl-^ - Rl'A {Rl_ +R\- Rl+ - Rl) 



2d 



Rt^R: 



R 



1+ 



Rf) {Rl_ + m - - R^, 



2d 



R{-JR 



(B2) 

(B3) 
(B4) 



which can be obtained from Eq. (Bl). Thus one can 
r estore the function R\2 in two steps: first solving Eqs. 



(B3 



B4 ) and then substituting the functions R\± into Eq. 



(Bl). Though at each step one should solve ordinary dif- 
ferential equations, this is a hard program and we are not 
able to perform it entirely. Nevertheless, we can examine 
some crucial cases. 

First of all, using the above scheme one can establish 
the behavior of the field R\2 where ri^2 are close to ±r/2. 
The analysis shows that the field R\2 is smooth near the 
points. This is the justification of our reduction proce- 
dure leading to the effective action (|3.(]|). 

Simple consideration shows that the most dangerous 
geometry t hat c ould lead to violation of the triangle in- 
equafities (2.15) is realized if two of the three points are 
close to ±r/2 whereas the third one (say, r\) lies in the 
middle between the points. In this case = R\- and 
the triangle inequality tells in this case that the difference 
2i?i_ — R^ should be positive. To prove the inequality 
let us write the equation for this quantity 



d(2i?i_ -RJ) 
di 

Dm^ [R 



D[2R\-J -R 



?l-7 



(B5) 



2-7 



R^S'')[R\ - R^_) {R--2R1- 



It is easy to show that at small times the quantity 
2Ri^ — R^ increases. Suppose that at some moment 
of time it becomes zero. Then, the second term in the 



R 



2-7 



right-hand side of (B5) is zero, while the first one is neg- 
ative. It means that the difference 2i?i_ — i?_ increases 
(recall that we move backwards in time) , though it should 
approach zero from the positive side. The contradiction 
proves that 2_Ri_ — i?_ is always positive. 

Let us now restore th e te rms subleading over \/d. 
Adding the terms to Eq. (Bl) we can repeat our consid- 
eration. Again, we should consider the same geometry 
as above. Writing the equation for 2i?i_ — we are 
convinc ed t hat in addition to the two terms presented 
in Eq. (B5) one should take into account also the term 
of the next order over \/d. This term is non-zero when 
2_Ri_ = therefore it starts to compete with the lead- 
ing term -DR'S^{2^ - 1) if 2i?i_ - i?_ is small. In this 
case the leading term is proportional to 2^ — 1 that is to 
7 at small 7. Therefore, if 7c? is not large, we can not 
make definite conclusion about t he sig n of the difference. 
Thus we arrive at the inequality ( 2.16 ) formulated in the 
main text. The crucial cases investigated above make us 
confident that the triangle inequality holds for our in- 
stanton in the whole space provided the inequality ( 2.16| ) 
is satisfied. 
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2. Fluctuations 



Here we extend the analysis of the fluctuations pre- 
sented in the Sec. 



II B that is suitable for all the points 



excluding vicinities of ±r/2. The reason is that the quan- 
tity V introduced by Eq. ( |3.11 ) is small during some 
stages of the evolution. Therefore, i?+ = R{r /2,r /2) 
and R- = R{r/2, —r/2) are almost equal to each other. 
Thus, we should check that fluctuations do not destroy 
this proximity. Having the problem in mind we will as- 
sume u <C 1 below. 



We can take the effective action (2.23) as the starting 
point of our analysis. It will be enough for our purpose 
to examine fluctuation effects in the harmonic approxi- 
matio n. Therefore, we should expand the effective action 
(2.23) up to the second order over the fluctuations SRJ2 
and 611112. The first order term of the expansion vanishes 
due to our saddle-point equations. The second order term 
can be written as 

il^^'^^i J dtdridr2Smi2l~^dtSR'l2 
D f 

-— / dt d'^r (3i2,345™i2<5m34 (B6) 
+ i / dtdridr2|y|VH^i2)/3(ri)/3(r2) 



J dtd'^r |25(3i2,34TOl2^TO34 + Qi2!34™12"^34| • 

where d^r = dri dr2 dr^ dr4. Here i?i2 is the field corre- 
sponding to the instanton (recall, that the way to restore 
it in the whole space was discussed in Appendix B 1 ) , 
mi2 is determined by the expression (3.2), and x^'^\ Q'"'^^ 
are the second order terms in the expansion over 6R^. 
As follows from Eqs. (3.1,3.3), the last two terms in (86) 
are relevant only for the points close to ±r/2. Therefore, 
in the general case w e can disregard these terms, and re- 
turn to the estimate (2.2(;), which can be obtained if only 
the two first terms in the action (B6) are kept. However, 



we are interested just in the behavior of i?i2 when the 
points Ti and r2 are close to ±r/2. In this case a special 
analysis is needed. 

First, note that short -scal e fluctuations of R12 are weak 
due to the restriction (2.15) which makes the amplitude 
of the fluctuations proportional to the scale. In other 
words, we should deal only with smooth functions R12. 
Next, due to the presence of the second term in the ef- 
fective action (B6) fluctuations of mi2 are relatively sup- 
pressed for points that arc not very close to ±r/2. That i s 
a consequence of the r-dependence of (3i2,34 [Eq. ( [2.11 )] 
which in the case w <C 1 has deep minima if the points 
are close to ±r/2 (more precisely some linear combina- 
tions have deep minima, they just determine the struc- 
ture of strong fluctuations of m). Therefore, relevant 
fluctuations of R12 can be estimated in terms of the ex- 



pression (3.3) and, consequently, in terms of the reduced 
action (3.6). Since in the main approximation over d 



the terms with ip2 and ips can be neglected in (3/7), the 
integration over to+ can be done explicitly. That leads 
simply to flxing th e exp ression (3.1C) and reduces the 
action to the form ( 3.13 ). 

So we should esti mate fluctuations of v and ji start- 
ing from the action ( 3.13 ) with the Hamiltonian ( 3.14 ). 
Actually, we should check the validity of the semiclas- 
sical approximation for this system with one degree of 
freedom. It is more convenient to perform this conven- 
tional procedure in terms of the canonically conjugated 
variables p and q wher e p = fJ.v/"f and v = exp(g). Then 
the Hamiltonian (3.14) is rewritten as 



H 



-IP- 



2(2-7) 



1 



7 



27 



exp((7) 



D7 



i?V^x'(i?+)exp((7), 



(B7) 



where we took into account w <C 1 . The subsequent anal- 
ysis shows that the semi-classical approximation is bro- 
ken only in the vicinity of the reverse point Vr that exists 
at n > ric- Let us estimate this vicinity. Near the re- 
verse point the last term in (B7) can be neglected and 
resolving the relation H = Hi we get 



1 ± 



7 



27 



where we substituted (3.26|). The semi-classical approxi- 



mation is broken Hp dp/dq 1 which gives 



Since the main contribution to the integrals such as Eq. 
( |3.2l|) is made at v ~ Vr ^ v,,, the above narrow vicin- 
ity (where the semi-classical approximation is broken) is 
irrelevant for our results. 

Of course, the above analysis of the fluctuations on the 
background of our instanton is not exhaustive. Never- 
theless, we believe that the arguments presented demon- 
strate the weakness of the fluctuations. 



APPENDIX C: 



As was mentioned in Sec. [11 A| , keeping only the main 
over 1/d term in the action (|3.6D is potentially dangerous 
because of a possible conflict of the limits d ^ 1 and 
L/r ^ 1. If this were the case, our results would not 
be applicable deep in the convective interval, where the 
ratio r/L is very small. Fortunately, this is not the case. 
In this appendix w e wU l consider the equations, follow- 
ing from the action ( |3.6| ), not neglecting terms subleading 
over 1/d. We show that r/L does not interfere with 1/d, 
and the scheme, presented in Sec. HI is reproduced with 
minor modiflcations. Repeating the scheme, we obtain 
the expressions for the anomalous exponents ^„ in this 
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formulation. To avoid a misunderstanding let us stress 
that our s chem e is correct only in the limit ^7 > 1 (see 
Appendix B 1 ) . Therefore, one could consider this ap- 
pendix only as a demonstration of the absence of the 
aforementioned conflict of limits. 

Below, we will assume that D = 1, L = 1, and P2 = 1 
which can be done by rescaling time t, the coordinates 
r, and the passive scalar 9. To restore the full answers 
one should add simply the factor P2/DL'' to i?^ in all 
expressions. 

Extrcmum conditions for the action ( |3.6D give four 
equations of motion for the quantities R± and m± . As 
before, the boundary conditions to the equations are 
i?+ — and R- = r at t = 0, and m± — > at t — > —00. 
The equa tions are canonical, with Hamiltonian given by 
Eq. (3.7). Since it does not depend on time explicitly, 
the energy is conserved. Moreover E — 0, which is a 
consequence of the asymptotic behavior m± of the 
instantonic solution at t —00. 

As it is known from classical mechanics, for a Hamil- 
tonian system which has an integral of motion one can 
reduce the number of degrees of freedom by one. In our 
case, we can pass from the system of four canonical equa- 
tions to that of two. Let us express say, m+ via the other 
variables with the help of the conservation law E ~ Q: 



m-^-Lp'i — a — m_(^2 



(CI) 



{m-tp2 - a)2 - ip3 [miifi + 2a{m- + \y\^U)\ , 



where 



u = x(R+) - x{R-) 



The sign in front of the square root in (CI) should be mi 



nus to ensure the correct behavior of i?+ at small time. 



Let us make the subst itution (3.11 ). Then we obtain 
a generalization of Eqs. ( 3.13 , 3T4^ 



H 



;x = y de Va«« - H] , 

a — fi(j)2 



-^(l + v) + 



(C3) 
(C4) 



- a)2 - ^3 + 2a [m + \y\^Uey^]} 

Note that 



We introduced here the notations 

(2(1 + i;)2/T- 1)21 



61 --(l + z;)2-4/7 
4(rf + 1 -7) 



1 



2-7 

62 = 2 • 



1 + v 

1 1 , 2+v 

^ 1 + v' 



(C5) 



(l + w)2/7-l (l + t;)2/7^ 



One can derive also a generalization of Eq. (3.12) 

no 

2na eT«C/dC 



y/{fi^2 - a)2 - h {/i^'/'i +2a[fi+ |2/|2[/e7?]} 



As before (see Sec. Ill A) we can divide the evolution 
into three stages. Let us analyze the first stage. Since 
U = there, the quantity H does not explicitly depend 
on ^ and therefore its value (which we designate Hi) is 
conserved during the first stage. Then from the relation 
we can express /i via v as follows 



H1S4 ih Fi 



Fi = jHfS-s + 2aHiS2 + a^v^ . 



We introduced the shorthand notations 



52 = 

53 = 
^4 = 



C^l + 2{1 + V)(j)2 + + 

01 +02 + (1 + V)(02 +03), 

02 - 0103 , 

02 + (l+w)03. 



Then we can derive an equation for v 



7 ^d^v — 



±SiFi 



XS2 + H1S3 T SiFi 



^ -Giv) 



(C6) 
(C7) 



(C8) 
(C9) 
(CIO) 
(Cll) 



(C12) 



This equation is the direct generalization of Eq. ( |3.19[ ) . 
To ensure the finite value of t he a ction (C3), one should 
t ake the lower signs in (|C^,|CT^). A solution of Eq. 
( C12 ) with the correct boundary condition is given by Eq. 
(OC), where G should be substituted from Eq. < \C12\) . 

Let us establish the behavior of G at small v which is 
important for description of the first stage. To do that, 
we should take into account the relations 



47-72 2 



-2-31- 



1)2 = 2 + v v% (j)3 = -2 + V - 

7 



^1 



Kd-7) 



^2 



8(d-7) 



16(d-7) 2 

63 = 5 : 



r 
2 



Si = — . 



valid at w ^ 1. Note that all the functions Si have the 
homogeneous behavior cx v'^ at small v. Therefore, start- 
i ng from r elation (3.20) we can find the same estimates 
( [3.24^^2^ ) for and Vr with small corrections of the 
order 1/d. 

Next, we should analyze the second stage and match 
it with the first stage. That is a repetition of the proce- 
dure described in Appendix |^. As a result we find the 
exponents 
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Cn = ^ (a - v/«' - 2a7n + 4(d - 7)^2) , (C13) 

(C14) 




In th e lim it 7^ ^ 1 we recover the previous results (3.3C) 
and (3.32). To avoid a misun derstanding, let us stress 
that the expressions ( C13| . C14 ) cannot b e used to estab- 
lish l/d corrections to the exponents ( 3.30| , 3.32) since 
contributions of the same order, related to the triangle 
inequalities ( 2.15| ) and fluctuations, are unknown. 



APPENDIX D: 

In this appendix we present a consideration of the sec- 
ond stage and matching conditions for the instanton so- 
lution. The procedure appears to depend strongly on the 
order n of the structure function. Therefore, we consider 
different cases separately . The designations used below 
were introduced in Sec. III. 



1. Intermediate tail 

Let us first consider the case where the behavior of the 
function is monotonic during the whole first stage. 
As was demonstrated in the main body of the text, this 
case is realized if Hi > He and that the variable v is 
small at the end of the first stage. We will show that 
the variable remains small also during the second stage, 
going to zero at the third stage. Therefore, the evolution 
of V during the last substage of the first stage and during 
the second and the third stages can be described in terms 
of the Hamiltonian (3.14), where the condition w ^ 1 is 
utilized, which simplifies the analysis es sentially. 

In the limit f <C 1 we get from Eqs. (3.14, |3.15 ) 



-A^- + ^^A^'-'-^<V(i?+)- (Dl) 



H = 



We see that the first two terms in the expression ( pi| ) 
depend only on fiv. This is the reason why the equation 
for the quantity 



\y\: 

D 



■<+V(i?+)^, 



(D2) 



that can be derived from Eqs. (Dl 3.16) contains in the 
right-hand side only the term proportional to x'iR+)- 
The term can be neglected if either i?+ ^ 1 or ^ 1 
that is during the first and during the third stages. 
Therefore, the quantity fiv is conserved there. Recall 
that due to the boundary condition /x = at ^ — > 00 and 
consequently fiv — during the third stage. 

Integrating relation (D2) from any ^ corresponding to 
the last substage of the first stage up to -l-oo we get the 
following integral relation 



D 



J d^Rl+^x'{R+)v, (D3) 



which is correct for t he fi rst stage if w ^ 1. It is instruc- 
tive to compare Eq. (D3) with the relation 



(D4) 



which can be obtained from Eq. (3.12 ) at u <C 1. Recall- 
ing also relation y = —in I'd (2.21), we get 



7n 



(D5) 



Substituting Eq. (D5) into (Dl) and neglecting the term 
with x' we find that the value of the Hamiltonian during 
the first stage is 



Hi 



in , 2-7 2 

n . 

2 2d 



(D6) 



During the second stage fiv diminishes from to 
zero. Therefore, the equation 



dint; / 4(2-7) 

7 -1+ , 2 > ^ 



^72 



shows that the quantity v does not vary essential ly during 
the second stage. Then we obtain from Eqs. (D4 2.21 ) 
the following estimates 



v^,L^ 



nP2 



\y\ 



L1P2V.^ 



(D7) 



Recall that is the value of v at the end of the first 

stage. 

Substituting Eqs. (|3.23| ) and into Eq. ( |3.1^ ) we 
get 



G = w (1 - n/ric) , 



(D8) 



where ric is defined by Eq. ( 3.2^ ) . Recall that the mono- 
tonic behavior of v implies G > and therefore the ex- 
pression in ( p8| ) is correct \i n < ric- It means that for 
n > Tic the instantonic solution of the considered type 
does not exist and we should look for another possibility. 
We postpone the problem up to the next subsection and 
continue to anal yze the monot onic v regarding n < ric. 
Substituting (p|) into Eq. (|3.24D we get 



In — = (1 — n/ric) 7ln — 
r 



(D9) 



Next, we should find the leading contribution to the ac- 
tion. This contribution is made mainly by the first stage 
produ cing a large logarithm. Therefore, we can write 
using (3.13) 



dv [ fj.+ 



Hi 
G 



16 



Substituting here Eqs. (D5,D(; .D8) we find with the same 
logarithmic accuracy 



2d r 



(DIO) 



Fi nally we can determine the s truc ture functions 5„ 
( L2| ) in accordance with formula ( ^.2[ ). Collecting Eqs. 



(D7.D9,D10) we obtain the expressions (3.29,3.30) of the 



main text. 



2. Remote tail 

Let us proceed to discussing the character of the in- 
stantonic solution at n > Uc- We can solve the corre- 
sponding equations in two limiting cases: n ^ ric and 
ji ~ Tie <^ ric. The latter case is considered in Appendix 
D 4 . In the subsection we accept ric, the inequality 
ensures also i;* 3> 1 (recall that v^, is the value of v at 
the end of the first stage). Due to the condition ^ 1, 
the last substage of the first stage and the second stage 
can be examined in terms of the Haniiltonian 



Hamiltonian (Dll) should be equal to Hi k, He, (3.22) 
that is 



At z)* ^ 1 the relation leads to 



/I* w 2c? ^ ^/^njd. 



(D15) 



Actually, /i « 2(i during the w hole last substage where 
3> 1. Deriving the expression (DIE) for w*, we used the 
estimates ( pl^ , pl4D . 

Now, we c an turn to the calculation of the effective 
s before, the main contribution to the 



action (3.13). As bef 
action is made by the first stage. Hence 



iX = y du + du 



(D16) 



Moreover, only the contribution related to the vicinity of 
the reverse point is relevant. Using relations (3.21) we 
obtain 



'fJ,V - 



2d' 



\y 



x{R+) - xiR+v'/') 



(Dll) 



which follows from Eqs. (|3.1l|j3.14| , |3.15|) at w > 1. At the 



end of the first stage R- ~ L whereas i?+ « R^/vV^ <C 
L. The second Eq. (3.16) shows that /i which is equal 
to zero at R- > L varies essentially at R- ^ L. To find 
the value of fi at R- < L, we will use the relation 



d 



[{H + nv)R'_^ 



71 _ 



\y\: 

D 



(D12) 



that follows from Eqs. ( |3.16| , |DT1| ). Actually ( pl2D is 
the equation for m+ that can be obtained from Eqs. 
( g!6|j3^ , |3!9| ) under the same conditions that led to Eq. 
( Dll JT^Since R+ <C R- ~ L during the last substage of 
the first stage and the second stage, the right-hand side 
of Eq. (D12) can be neglected and we get a conservation 



law for {H + fiv)R^^ . Equating the values of that quan- 



tity at _R_ < L and at R- 
of the first stage 



> L we find that at the end 



d\y\^P2 



D 



(D13) 



Now, let us return to relation (p.4|). Since i?_ ^ L ^ 
at the end of the first stage the main contribution to 
the integral (3.4) is made when already R- > L and R+ 
increases to L. Using Eqs. ( ^^4 3.1C ), and ( 2.21 ) we get 



(D14) 



Substituting the expression (D14) for \y\ into (D13) we 



conclude that with n. In addition, at i?_ < L the 



r 



(D17) 



Here, we substituted Hi « He and took into account that 
/i has no sin gular denominator which is clear from the re- 
lation ( |3.17| ). And finally, we can determine the structure 
functions 5„ ( |1.2D in accordance with (|3.2| ). Collecting 
( D14 D17) we obtain the expression ( p.33D of the main 

text. 

Expression (|D15|) shows that indeed ^ 1 (which is 
the applicability condition of the above consideration) at 
n 3> if there are no additional small parameters. The 
next subsection is devoted to special situations appearing 
in the presence of such parameters. 



3. Special cases 

In the subsection we treat three special cases that are 
realized at n > ric if 7 ^ 1 , 2 — 7 ^ 1 , or n is e xtre mely 
large. The general picture given in Appendix D 2 does 
not essentially change, but particular answers should be 
slightly corrected. 

The case 7 <C 1 needs a special consideration since as 
follows from Eq. ( D15 ), the condition ^ 1 is violated 
at n ~ (i/7 3> nc ~ d'y. Therefore, an intermediate re- 
gion exists at Uc - C n <C (i/7 where 7 <C <C 1. Then 
we get from Eq. ( 3.14 ) 



H 



where 



EL 

2d 



)~x{R-)], (DI8) 



R- = R+ exp (w/7) 



(D19) 
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Equations ( [3.16 ) now read 



dv /i 



■RlR^x'iR-) 



(D20) 
(D21) 



At the e nd of the first stage i?_ ^ L. As follows from 
relation ( D19 ), due to ^ 7 the estimate i?+ <C i is 
valid. Therefore, the integral (3.4) for is determined 
by the time interval between the moments when R-^- and 
i?_ go through L. We can aga in pa ss to integrating over 
i?+ in accordance with Eq. ( 3.10| ). Though i?+ ^ L 



at R- ~ L, 'jR'l « i?l(l — ti) due to smallness of 7. 
Therefore, we return to the estimates (D7). Next we 
can estimate /i* from Eq. (D21) where the term 7/i is 
dropped. Since i?^ ^ L'' at R_ ~ L, integrating over ^ 
we obtain from Eq. (D21) 



D 



(D22) 



To justify the estimate ( P22| ), one should check that v 
does not va ry es sentially in the region, where i?_ ~ L. 
Using Eq. ( D20 ) we obtain that the c ondit ion is satis- 
fied. Substituting formulas (^) into ( p2^ ) we obtain 
the previous relation fi^v^ ^ 3> \Hc\. Equating then 
the Hamiltonian (D18) to He at v — we get 




(D23) 



The estimate ( |D23| ) for s urpr isingly coincides with 
( pT^ ). Next, we obtain from ( p7| ) 



71c D 



(D24) 



Therefore, the expression (|3.33| ) should be corrected and 
we get 



n P2C3 



rir D 



V2 



(D25) 



where C3 is a non-universal parameter of order unity and 
Cc is defined by Eq. (|3|). 

Some peculiarit ies ap pear also at 2 — 7 ^ 1. In this 
case the function ( 3.15D can be approximated as 



(j) « 2dv In V , 

i f 7j ^ 1 but (2 — 7) In w ^ 1. Then we obtain from Eq. 



H 



-fiv 
D ^+ 



2-7 
2d 



(D26) 



Expression (D26) has the same structure as (Dll) except 
for the logarithmic factor. Taking into account u 3> 1 we 
conclude that the Hamiltonian (D26) leads to the con- 
servation law ( |D12| ) which is satisfied with the same ac- 
curacy l/v as previously. Therefore, we get instead of 



d\y\^P2 



D{2--i)\nv. 



nd 



(2-7) ln^;* 



(D27) 



where we used ( D14 ), its validity is accounted for the 
inequality v ^ 1. O ne can check that the main contribu- 
tion to H in (D26) at the end of the first stage is made 
by ~fiv. Equating —fi^v^ to He we get 

TL 

Inw, . (D28) 

71c 

Now we obtain the condition 71 <^ 7ic/(2 — 7), under 
which the regime under consideration is realized. Since 



V ^ 1 in the regime we have the same expr essio n (D14) 
for and, consequently, the same estimate ( 3.33 ) for the 
structure functions. 

We finish this subsection with a discussion concern- 
i ng ex tremely large n. In accordance with the estimates 
(D15) Va. increases with increasing n. Therefore, the 
last substage of the first stage where v ^ 1 starts to 
play a n essential role. Fir st of all we should correct Eq. 
( |3.2lD returning to Eq. (HH). U v = there then 
in the right-hand side of the relation can be sub- 
stituted by . The contribution to the integral in the 
left-hand side of Eq. (3.21) associated with the last sub- 



stage (when V increases from 1 to w*) can be found sub- 
stituting G{x) « ~x there. Thus the substage makes the 
additional logarithmic contribution to the integral. Tak- 
i ng th e contribution into account we get instead of Eq. 

(P3) 



27r 



8(2-7) ^{He-H,) 



7 In ■ 



L 



21nt;* 



(D29) 



There appears also an analogou s cor rection to expression 
( D17 ). Calculating the action ( |D16| ) we get 



iJ w Hr In 



L 



lnu« 



(D30) 



Here we neglected the term originated from J fidv since 
using (D15) the term can be estimated as \/n/ (7c?) and 
is consequently negligible in comparison with the term s 
^ n e nter ing Sn via Substituting the expression (D3C) 
into (3^) we conclud e tha t the correction related to the 
last contribution in (D3C) can be neglected in compar- 
ison to the strong 71-dependence of 7?". Nevertheless, 
the presence of th e additional term in the right-hand 
side of Eq. ( D29| ) shows that our approach is broken 
at 21n7;* ~ 7ln(_L/r) since He — H i ceases to be a small 
there. Using the estimates (D15) we arrive at the con- 



dition (3.35), which is the applicability condition of our 
approach. Note that at \n(n/d) ^ 7ln(L/r) a scaling 
behavior of Sn is destroyed. 
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4. Instanton near critical order 

Here we wil l con sider the case when n is close to the 
critical value ( [3.28 ): |n — nd <C jt-c. Then <C 1 and at 
examining the last substage of the first stage and of the 
second stage we can use the Hamiltonian, expanded over 
small V. 

li n < Tic then we can take the Hamiltonian (Dl). The 



only peculiarity of the consideration near is th at the 
logarithmic contributions to the effecti ve ac tion (3.12) 
and to the left-hand side of the integral ( 3.20 ) depend on 
n — Uc- Indeed, the expression (3.25) shows that G{x) is 
a linear function of x only in the restricted interval 



■ 7 



2d7 



He - Hi \ > V > 



Just the interval determines the value of the logarithm. 
I t means that we should substitute in the expressions 



In- 



In 



{n - UcY 



The substitution does not change the final answer (DlC) 
for the action but influences the value of d"^ because of 
( D7|). Taking the fa ct int o account we get the expression 
( P^6D that replace ( |3.29D . 

It is clear that at n Uc the difference n — ric starts to 
compete with r/L what destroys our construction. Let 
us estimate the corresponding value oi ric ~ n regarding 
n < ric- It follows from the asymptote ( |3.25| ) that our 
consideration is correct if {He — S> d'yv^/{2 — 7) . 
Expressing in the inequa lity v * via r/L from Eq. (3.24) 



we obtain the criterion ( |3.37 ) written in the main text 
Note that the criterion implies the inequality 

, L 
7ln - > 1 , 
r 

which should be assumed for our treatment to be correct. 
Otherwise a special consideration is needed. 

Let us now proceed to the case n > ric- We should 
take into account the next term of the expansion of H 



over V in c omparison with (Dl) which should correct the 
value ( p^ ) of Hi since we know that Hi is equal to H^,- 
The corresponding expression is 



H 



2(2-7) 
^72 



1 '^^^ \ 2 2 
1 — V 11 V 



27 



[x{R+)-x{R-)] 



(D31) 



where the difference x(i?+) — x{R~) shou ld be expanded 
up to the second order over v. Eqs. ( 3.16 ) with the 
Hamiltonian (D31) lead to 



lAv 2, 
--rz = -V + fj,v [ai 

7 



1 d 



a2v) , 



(D32) 



7 dC"^^"^ " + ^i?X+'x'(i?-)« , (D33) 

2(4-7)(2-7) tT^',A^ 
"2 = j:^, ■ (D34) 



ai 



4(2-7) 
^72 



Equating the Hamiltonian (D31) to the value ( 3.22| ) we 
get 



2 



(D35) 



at the end of the first stage. Note that due to <C 1 we 
can use the estimates (D7). 

Our next aim is to estimate in terms oi n — ric- To 
do it, we will use the identity dH/d^ — dH/d^ that is 
correct for any canonical system. Let us integrate the 
relation over the second stage. Since H — Hi at the first 
stage and H = at the end of the second stage we obtain 

-1/1 = ^1 de {7^ [x(i?+) - x{R-)] 

+Rl[R+x'{R+)-R-x'iR-)]}- (D36) 

Expanding the right-hand side of relation ( |D36| ) up to the 
second order over v, integrating b y parts (to remove high 
derivatives of x) and using Eqs. ( 3.12 . D32 , D33 , 2.21 ) one 
can obtain the expression 



'Hi = — [I — j - — (n-nj— / dS^R^^x iR+)f^v' 

/ d^Rl^'x'iR+ W-^aia, I d^Rl+''x'iR+)v / d^^'v' 



(D37) 



The first term in the right-hand side of (D37) reproduces 
the previous result (D6) and other terms are small cor- 
rections proportional to u*. 

We know that in the main approximation over r/L we 
can substitute Hi = He- Regarding n — Uc <C nc we ne- 



glect the second term in the right-hand side of Eq. ( D37 ) . 



The last two terms in the right-hand side of Eq. (D37) 
can be estimate d us ing ~ jric, x ^ P2, R+ ~ L and 
the estimation (D7) for jyp. Combining all together we 
get 



(D38) 
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Let us stress that d ue to a2 > the last two terms in the 
right-hand side of ( D37 ) are positive if x(i?) is a mono- 
tonic function which is a reasonable condition. Therefore, 



a solu tion of Eq. ( D37 ) for w*, determined by the esti- 
mate (D38), really exists. We conclude that there is an 
instantonic solution with non-monotonic behavior of v at 
any n — Uc <^ fic- 

The last assertion should be corrected since it is true 
only if L/r — > cx). At finite L/r there is a narrow region 
of very small n — ric where our scheme does not work. To 
establish the corresponding criterion let us remind that 
our consideration is valid if Vr <C . It follows from 
Eq. ( |3.27| ) that ~ 1 / [-y (L/r)] . Then, from formula 
( D3^ ) we get the inequality ( 3.37 ). 
And fi nally , we can determine the structure func- 



tions Sn (|l.2 ) in accordance with (3.2 ). Co llecting Eqs. 
(D7.D17,D3S) we obtain the formula ( 3.36 ) of the main 
text. 



[1 

[2: 
[3; 
[4: 
[s: 
[6; 

[r 

[8 

[9: 

[lo; 

[11 

[12 
[13 

[14 



A. N. Kolmogorov, C. R. Acad. Sci. URSS, 30, 301 
(1941). 

H. L. Grant, R. W. Stewart, and A. Moilliet, J. Fluid 
Mech. 12, 241 (1961). 

G. K. Batchelor and A. A. Townsend, Proc. R. Soc. Lon- 
don Sect. A 199, 82 (1962). 

C. Meneveau and K. R. Sreenivasan, J. Fluid Mech. 224, 
429 (1991). 

G. K. Batchelor, Theory of Homogeneous Turbulence 
(Cambridge University Press, New York, 1953). 
A. Monin and A. Yaglom, Statistical Fluid Mechanics 
(MIT Press, Cambridge, 1975). 

U. Frisch, Turbulence: the Legacy of A. N. Kolmogorov 
(Cambridge University Press, New York, 1995). 
A. M. Obukhov, C. R. Acad. Sci. URSS, Geogr. Geo- 
phys., 13, 58 (1949). 

S. Corrsin, J. Appl. Phys., 22, 469, (1951). 

R. Antonia, E. Hopfinger, Y. Gagne, and F. Anselmet, 

Phys. Rev. A 30, 2704 (1984). 

C. Meneveau and K. R. Sreenivasan, Phys. Rev. A 41, 
2246 (1990). 

K. R. Sreenivasan, Proc. R. Soc. (London) A 434, 165 
(1991). 

L Hosokawa and K. Yamamoto, in Turbulence and Coher- 
ent Structures, ed. by O. Metais and M. Lesieur (Kluwer, 
London 1991). 

R. H. Kraichnan, Phys. Fluids 11, 945 (1968). 



[15 
[16 
[17 

[is: 

[19 

[20 

[21 

[22; 

[23; 

[24; 

[25; 

[26 
[27 
[28' 
[29 

[30 
[31 
[32; 

[33 

[34' 
[35; 

[36 

[37 

[38' 
[39 

[40 
[41 

[42; 

[43; 

[44; 



B. Shraiman and E. Siggia, C. R. Acad. Sc. 321, Ser. 11, 
279 (1995). 

K. Gawedzki and A. Kupiainen, Phys. Rev. Lett. 75, 
3608 (1995). 

M. Chertkov, G. Falkovich, L Kolokolov and V. Lebedev, 
Phys. Rev. E 52, 4924 (1995). 

M. Chertkov and G. Falkovich, Phys. Rev. Lett. 76, 2706 

(1995) . 

D. Bernard, K. Gawedzki and A. Kupiainen, Phys. Rev. 
E 54, 2564 (1996). 

B. Shraiman and E. Siggia, Phys. Rev. Lett. 77, 2463 

(1996) . 

E. Balkovsky and D. Gutman, Phys. Rev. E 54, 4435 

(1996) . 

E. Balkovsky, G. Falkovich, and V. Lebedev, Phys. Rev. 
E 55, R4881, (1997). 

R. H. Kraichnan, V. Yakhot, and S. Chen, Phys. Rev. 
Lett. 75, 240 (1995). 

A. L. Fairhall, B. Galanti, V. S. L'vov, and L Procaccia, 
Phys. Rev. Lett. 79, 4166, (1997). 

U. Frisch, A. Mazzino, and M. Vergassola, submitted to 
Phys. Rev. Lett. (1998). 

R. H. Kraichnan, Phys. Rev. Lett. 52, 1016 (1994). 

V. Yakhot, Phys. Rev. E 55, 329 (1997). 

M. Chertkov, Phys. Rev. E 55, 2722 (1997). 

P. C. Martin, E. Siggia, and H. Rose, Phys. Rev. A8, 423 

(1973). 

C. de Dominicis, J. Physique (Paris) 37, cOl-247 (1976). 
H. Janssen, Z. Phys. B 23, 377 (1976). 

E. Balkovsky and V. Lebedev, submitted to Phys. Rev. 
Lett. 

G. Falkovich, L Kolokolov, V. Lebedev and A. Migdal, 
Phys. Rev. E54, 4896 (1996). 

V. Gurarie and A. Migdal, Phys. Rev. E54, 4908 (1996). 
E. Balkovsky, G. Falkovich, I. Kolokolov and V. Lebedev, 
Phys. Rev. Lett. 78, 1452 (1997). 

G. Falkovich and V. Lebedev, Phys. Rev. Lett. 79, 4159 

(1997) . 

E. Balkovsky and G. Falkovich, Phys. Rev. E 57, 1231 

(1998) . 

B. Shraiman and E. Siggia, Phys. Rev. E 49, 2912 (1994). 
L. Richardson, Proc. Roy. Soc. (London) A 110, 709 
(1926). 

A. N. Kolmogorov, J. Fluid Mech., 12, 82 (1962). 
M. Lesieur, Turbulence in Fluids (Kluwer Academic Pub- 
lishers, Dordrecht, 1997). 

D. Bernard, K. Gawedzki and A. Kupiainen, J. Stat. 
Phys. 90, 519 (1998). 

U. Frisch, A. Mazzino, and M. Vergassola, in Proceedings 

of Europ. Geophys. Soc. meeting, 1998. 

R. H. Kraichnan, Phys. Fluids, 9, 1937 (1966). 



20 



